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Abstract 

We study the single transverse spin asymmetry for the inclusive pion production 
in the nucleon-nucleon collision, p^p — > irX, based on the twist-3 mechanism in the 
collinear factorization. We derive the soft-fermion-pole (SFP) contribution to the 
twist-3 single-spin-dependent cross section associated with the twist-3 quark-gluon 
correlation functions in the polarized nucleon. We find that the SFP can give rise to 
a large effect to the asymmetry owing to the large partonic hard cross sections 
with the large color factor, if the SFP function has a similar magnitude as the soft- 
gluon-pole (SGP) function, in spite of the absence of the "derivative term" for the 
SFP function unlike for the SGP function. 
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The large single transverse spin asymmetry (SSA) observed in p^p —>■ irX at FNAL 
(V^ = 20 GeV) [1] triggered lots of theoretical and experimental studies to clarify the 
origin of the asymmetry. (See [2] for a review.) In the conventional framework for the 
high-energy inclusive reactions, i.e., parton model and perturbative QCD, the asymmetry 
An = (cr^ — 0"^)/ {a^ + a^) is of the order of agmg/px where pt is the transverse momentum 
of the pion and niq is the mass of u and d quarks, and hence had been expected to be 
negligible [3]. The transverse polarization of hyperons produced in unpolarized pp collisions 
was another example of the unexpected large SSAs in early days [1]. More recently the large 
An was also found at higher energy at RHIC (v^ = 200, 62.4 GeV) [SlElEj. HERMES and 
COMPASS also reported large SSA in semi-inclusive deep- inelastic scattering (SIDIS) [SlIS]- 

Understanding of these SSAs requires the extention of the framework for QCD hard 
processes, and by now the mechanisms leading to such large SSA have been understood 
to some extent in terms of these new frameworks, which are often classified into two cat- 
egories. One is the so-called "T-odd" distribution and fragmentation functions [TOl [TT] in 
the transverse momentum dependent (TMD) factorization approach, and the other is the 
twist-3 quark-gluon correlation functions in the coUinear factorization approach p!3l (Til fT5] . 
The former approach has had phenomenological successes in the description of SSAs in var- 
ious processes [12]. Universality property of the TMD functions have been also studied in 
detail [Ml [171 [HI [iHl [20], while the factorization proof has been limited to e^e~, SIDIS 
and Drell-Pan processes [211 [22l [23] . In the framework of collinear factorization SSA is a 
twist-3 observable, and this approach is applicable to wider classes of processes including 
p^p — > nX, although its validity is limited to large pt region such that px can be regarded 
as "hard". The formalism of calculating twist-3 single spin dependent cross section with 
twist-3 distributions has been developed in [TH |T5], in particular, gauge invariance and 
factorization property of the cross section was proved in [15] in leading order QCD. So 
far the formalism has been applied to SSA in Drell-Yan [H [Ml [27], SIDIS [151 [25l [26l [27] , 
^ ttX [281 [291 [301 [an [32], -> A^X [331 [Si [35] , heavy-quark production in pp colli- 
sion [361 [2Z] etc. Since SSA is a naively "T-odd" observable, it appears as an interference 
between the scattering amplitudes which have different complex phases. In the twist-3 
mechanism of SSA, this phase arises as a pole part of a relevant internal propagator in 
the hard part, and those poles are classified as soft-gluon-pole (SGP), soft-fermion-pole 
(SEP) and hard-pole (HP). It is also shown in [271 [22] that the SGP contribution can be 
related to a certain twist-2 cross section in general. The connection between the TMD 
approach and the twist-3 approach at the intermediate region of the transverse momentum 
was studied, and it's been shown that these two mechanisms provide a unique and con- 
sistent QCD description of SSA for the Drell-Yan and SIDIS [Ml [2S1 [M]- Based on these 
developments, we will study in this letter the SSA for p^p —>■ nX in the framework of the 
collinear factorization. 

The cross section for pp — ^ nX can be written as a convolution of two distribution 
functions associated with the initial protons, fragmentation function for the final pion, and 
the corresponding partonic hard cross sections. In the twist-3 cross section, one of three 
nonperturbative functions becomes a relevant twist-3 function, and thus the single-spin- 
dependent cross section contributing to p^p nX consists of three kinds of twist-3 cross 
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sections: 

Af7*"3 ~ G^^\xi,X2)0f{x')®D{z)®aA 

+ 5q{x)® f{x')®&^\zi,Z2)®ac, (1) 

where the first, second and third factors in each contribution are, respectively, the polarized 
distribution function in p\ unpolarized distribution in p and the fragmentation function 
for the pion, and cfa,b,c are the partonic hard cross sections. The two-variable functions 
G'^^\xi,X2), E^'^\x'i^x'2)-, E'^'^\zi^Z2) are the twist-3 functions, and 5q{x) is the transversity 
distribution in the nucleon. f\x') and D{z) are, respectively, usual twist-2 unpolarized 
(quark or gluon) distribution and fragmentation functions. The analysis of the first, second 
and the third term of ([T]) has been performed in Refs. [281 El] , [29] and [30] , respectively. In 
these analyses, only the SGP contributions have been included. Kouvaris et al. derived the 
complete SGP cross section for the first term and presented a phenomenological analysis 
of FNAL and RHIC data [31]. The second contribution was shown to be negligible due to 
the smallness of the hard cross section ctb [21] • This is because only particular diagrams 
are allowed to contribute due to the chiral-odd nature of the transversity and E^'^\ which 
remains unchanged by the inclusion of the SFP contribution. In [30] , the so-called "deriva- 
tive" term of the SGP contribution from E^^\z^z) in the third contribution was analyzed, 
and it was shown that this contribution can be as large as the first one if the nonperturba- 
tive function E^^\z^z) has a similar order of magnitude as G^'^\x^x). Recently, however, 
it's been shown [^ that the SGP function E^'^\z^z) becomes identically zero, and thus 
the third term in ([1]) requires reanalysis including the non-partonic pole contribution from 
(i.e. imaginary part of) E^^\zi^Z2)- In this circumstance, we shall focuss in this letter on 
the analysis of the first contribution in ([1]), in particular, we will derive the complete cross 
section formula including the SFP contribution together with the SGP cross section derived 
in [21]. 

There are two independent twist-3 quark-gluon correlation functions (for each quark- 
flavor a), G'^{xi,X2) and G^(a:i,a;2), in the transversely polarized nucleon contributing to 
the first term in ([TJ.They are defined as (see eg. [T5] ) 

M^^{x,,X2) = 1^1 ^e^'^^e^'^(^-^^)(P5x|^J(0)^7F°^(/xn)n^C(Ari)|P5x) 

= ^Wi^e^'^^'^GUx,, X2) + ^^m)^,SlG''F{^^, x^) + . . . , (2) 

where is the nucleon mass, n is a lightlike vector satisfying p ■ n = 1, e"P'"'^-L = 
^yxP^n" S]_ with 60123 = 1 and + . . . denotes twist-4 or higher. By T- and P-invariance, 
two functions have the symmetry property under Xi ^ X2 as 

G^(xi,X2) = G^(x2,xi), G^(xi,X2) = -G^(a:2,xi). (3) 

With the gluon's field strength gF'^^np replaced by the covariant derivative -D" = d°' — 
igA°' in ([2]), one can define another set of twist-3 distributions. The relation between 
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those functions and the above {Gp,Gp} has been clarified in |25]. The twist-3 correlation 
functions for the "anti-quark" flavor and Gp can be defined from ([2]) by replacing the 
nonlocal operator xp'^ (0) g F°'f^ {nn)npijj"' (Xn) by its charge conjugatioE0, and is related to the 
original quark-gluon correlation functions in as 

GUxuX2) = G^(-X2,-xi), GUxi,X2) = -GU-X2,-xi). (4) 

For p^p nX, SSA occurs from SGP and SFP contributions, which, respectively, set 
Xi = X2 and Xj = (i = 1 or 2). Due to the above symmetry property, Gp contributes 
through SGP and SFP, while Gp contributes only through SFP. Thus the general structure 
of the twist-3 cross section takes the following forn|§: 

+ GFiO,x)0 f{x')0D{z)0asFP + GFiO,x)0 f{x')0D{z)0asFp, (5) 

where we have used the fact that the SGP contribution appears in the combination of 
Gf{x, x) — x^^^^^ [31], [32]. The origin of this combination and the connection of (Jsgp to 
the twist-2 cross section was clarified in [32]. Previous analyses focussed on the SGP term 
assuming it to be a dominant contribution, in particular, based on the observation that the 
SGP function receives enhancement by the "derivative" at large x [2H]. However, there is no 
clue that the SFP functions Gf{0, x) and Gf(0, x) themselves are small, and its importance 
depends on the behavior of the partonic hard cross sections (Tsfp and (5"sfp, even though it 
does not appear with the "derivative" . The purpose of this letter is to investigate the SFP 
contribution to p^p — > nX. We will derive a complete twist-3 cross section formula for the 
SFP contribution and study its impact on A^- in comparison with the SGP contribution, 
assuming that the SFP functions Gf{0,x) and Gf{0,x) have the same order of magnitude 
as the SGP function Gf{x,x) [iU] . 

We first recall the formula of calculating twist-3 cross section for p\p,S±) + p{p') — »■ 
Tx{Ph) +X following [15]. The contribution to the first term in ([1]) arises from the diagrams 
shown in Fig. 1, where the distribution function for the unpolarized proton (bottom blob) 
and the fragmentation function (middle small blob) for the pion are already factorized. Thus 
the parton lines coming from the unpolarized proton have the momentum x'p' and the one 
fragmenting into the pion have the momentum Ph/z, which are collinear to the momenta 
of the parent hadrons. The quark and gluon lines with the momenta /ci, /c2 and /c2 — ki are 
those from the polarized nucleon (top blob). The partonic hard parts Sp{ki,k2,x'p',Ph/z) 
in Fig. 1 (middle large oval) are defined for the nucleon matrix element ~ {pS ±\iIjA''iIj\pS ±) . 
As was demonstrated in [T^], in order to get the gauge invariant twist-3 cross section for 
SSA, one needs to reorganize the collinear expansion for Sp{ki,k2,x'p',Ph/z) with respect 

^Here we adopt the convention that the gluon's field strength gF"'^ is also replaced by its charge 
conjugation (together with ■0^-0'^ = C'ijr ) when we define the correlation fmiction for the "anti-quark" 
flavor. 

^Here and below we often suppress the flavor indices from the distribution and fragmentation functions. 
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Figure 1: General structure of the cut diagrams for the twist-3 cross section for p'p irX 
associated with the twist-3 quark-gluon correlation function in the polarized nucleon. Left 
(right) diagram represents contribution with the unpolarized quark (gluon) distribution in 
the unpolarized nucleon. 



to ki and /c2, and one eventually obtains the twist-3 cross section as 



dx' f dz 



X 



dxi / dx2Ti 



iu%M^{xi,X2 



dSp{h,k2,x'p',Ph/z)pP 



xf{x')D{z), 



(6) 



where Tr[---] denotes trace over spinor and color indice^, Mp(xi, X2) is the correlation 
function defined in ([2]) and u;"^ = g'^p — p'^rtp. All the other terms which contributes in the 
twist-3 level cancel out due to the following two relations [15] 



(X2 - Xi) 



dS,ik,,k2,x'p',Ph/z)pP 

dk^ 



-Sa{XiP,X2P,x'p',Ph/z) = 0, 



(7) 



and 



dSp{ki,k2,x'p',Ph/z)pP 

dkf 



hj_ — x-ip 



dSp{ki,k2,x'p',Ph/z)pP 

dk? 



0. 



(8) 



hi — Xip 



As noted before, single-spin-dependent cross section occurs from the pole (SGP and SFP) 



contribution in 



agp(fci,fc2,xV,Pfc/z)pp 



dkS 



hi — XiP 



The pole parts of the sum of the diagrams giving 



rise to SGP and SFP separately satisfy these relations, in particular, ([7]) is a direct conse- 
quence of the Ward identity for the pole parts of the diagrams (see [151 IH] for the detail). 



Color indices connecting Mp and Sp in ^ are implicit. 
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(a) 




(b) 

Figure 2: Typical examples of the cut diagrams for the SFP hard part, S^^^ in ([9]). In (a), 
two quark- lines coming from the polarized nucleon are in the opposite side of the cut, while 
in (b) they are in the same side of the cut. The bared internal line in each diagram generates 
SFP. Momentum for each line is defined as in Fig. 1 with ki = Xip. In this example, the 
coherent gluon line can attach at three different places of each diagram without it, and 
those three are represented by a single one shown in the right where the positions of the 
gluon-attachment are represented by the numbered dots (1),(2) and (3). 



For the case of the SFP contribution, which is our main interest in this paper, xi ^ x^ and 
one can use ([7]) in ([2]) to obtain 



Eh 



SFP 



API 



S 



xTr 



dx' 
x' 



dz 



dxi / dx 



Xi - X2 



iuj%M^{xi, X2)S^^''' (xip, X2P, x'p', Ph/z) 



^SFP 



f{x')D{z) 



(9) 



where {xiP,X2PtX'p' , Ph/ z) is the hard part for the SFP contribution. This formula 
guarantees absence of the "derivative" term for the SFP contribution and provides a simpler 
method for the calculation compared with iQ. 

To obtain the SFP cross section, one needs to calculate the diagrams with an extra 
coherent-gluon line as shown in Fig. 1 with ki = Xip {i = 1, 2). Those diagrams are 
classified into two types, typical examples of which are shown in Fig. 2. Fig. 2(a) shows 
the diagrams where the extra gluon line is attached to the cut diagram representing a twist- 
2 hard cross section. The line with a bar denotes the propagator generating a SFP. For 
each twist-2 diagram, there are three ways of attaching the coherent gluon, and thus we 
denote the sum of those three diagrams by a single diagram as shown in the right side of 
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Figure 3: Diagrams for the SFP contribution of the type shown in Fig. 2(a) in the quark (or 
anti-quark) fragmentation channel with the unpolarized quark (or anti-quark) distribution. 
The propagator with a bar gives an SFP. The parton c fragments into the final pion. The 
solid vertical line without an arrow represents the final state cut. The mirror diagrams also 
contribute. 

Fig. 2(a) where three dots on the internal lines represent the places to which the gluon line 
is attached. The mirror diagrams of Fig. 2(a) also give rise to the SFP contribution, and 
should be included in the calculation. Each of those mirror diagrams give the same result 
for the SFP cross section as the original diagrams. 

The diagrams shown in Fig. 2(b) have the quark and anti-quark lines from the polarized 
proton in the same side of the final-state cut, and the coherent gluon line is attached to the 
other side of the cut. Most of this type of diagrams are obtained by shifting the position of 
the final-state cut in the diagrams shown in Fig. 2(a). For these diagrams, there are again 
three ways of attaching the coherent gluon line, and we denote the sum of those three by 
a single diagram with three dots shown in the right side of Fig. 2(b). For the SFP cross 
section, some contributions in the same channel from the two types of diagrams shown 
in Figs. 2(a) and (b) cancel each other, as was the case for SIDIS [3B], which makes the 
calculation simpler |41j . 

Figs. 3~8 show the diagrams for the SFP contributions {S^^ {xip,X2P,x'p' , Ph/ z) in 
([9])). Three or four dots numbered as (1) ~ (4) in each diagram represent places where the 
coherent gluon is attached as explained in Fig. 2. Some diagrams in the figures have 4 dots, 
since the coupling of the coherent gluon through the 4-gluon vertex has to be taken into 
account for those diagrams. In Figs. 3~8, we have only shown the diagrams which represent 
SFP partonic hard cross sections for the "quark" distribution G°p{Q,x) and G^(0,a;) with 
X > 0. By reversing the direction of the arrows of all the fermion lines in each diagram, 
one obtain the hard cross sections for the "anti-quark" twist-3 distributions defined in (jl]). 

Figs. 3~6 represent diagrams in the channels with quark or anti-quark distribution for 
the unpolarized nucleon. Fig. 3 represents diagrams of the type Fig. 2(a) and Fig. 4 
represents those of the type Fig. 2(b) in the quark (or anti-quark) fragmentation channel. 
Fig. 5 represents diagrams of the type Fig. 2(b) and Fig. 6 represents those of the type 



7 



(a) (b) (<:) 




Figure 4: Diagrams for the SFP contribution of the type shown in Fig. 2(b) in the quark (or 
anti-quark) fragmentation channel with the unpolarized quark (or anti-quark) distribution. 
See the caption to Fig. 3. 




(f) (9) 



Figure 5: Diagrams for the SFP contribution of the type shown in Fig. 2(b) in the gluon 
fragmentation channel with the unpolarized quark (or anti-quark) distribution. See the 
caption to Fig. 3. 
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(2) 



"0W0"0"0" 

(3) 

(a) 







(f) 



Figure 6: Diagrams for the SFP contribution of the type shown in Fig. 2(a) in the gluon 
fragmentation channel with the unpolarized quark (or anti-quark) distribution. See the 
caption to Fig. 3. 

Fig. 2(a) in the gluon fragmentation channel. Fig. 7 represents diagrams with the gluon 
distribution for the unpolarized nucleon and quark (or anti-quark) fragmentation function 
for the pion. Fig. 8 represents those with the gluon distribution for the unpolarized 
nucleon and the gluon fragmentation function for the pion. In Figs. 3~8, we have omitted 
the diagrams which cancel each other as was observed for SIDIS [38l HI] . 

By calculating all these diagrams, we have eventually obtained the following cross section 
for the SFP contribution 



Eh 



^s^^SFP ^ ^2 ^pnP^Sj_ dz dx' f dx 1 



S 2 



X' 



X x'S + T/z 



6\ X 



-x'U/z 
x'S + T/z 



J2 (C^O, x) + G^(0, x)) {q\x') {D%z)aab^, + D\z)aai,^,) 

.a,b,c 

+q\x') {D\z)a,i_, + D\z)a,i_,)] 

a,b 

+ (g'f(0, x) + G^(0, x)) G{x') {D' 

a,c 



[z)aag^c + D^{z)a 



(10) 



+ (g^f(0, x) + G^(0, x)) G{x')D\z)aag-. 

a 

where the Mandelstam variables for the process are defined as 5 = {p + p'Y = 2p ■ p', 
T = {p- Phf = -2p ■ Ph and U = {p' - Phf = -2p' ■ Ph, q\x') denotes unpolarized 
quark distribution for the quark flavor 6, G{x') denotes unpolarized gluon distribution, and 
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D^{z) {i = c for the quark-flavor and i = g for gluon) is the pion fragmentation function. 
In the summations ^ ^, ^, ^ and the sum of a should be taken over all 
quark and anti-quark flavors (a = u, d, s,u,d,s,---) and the sum of b and c is restricted 
onto quark flavors when a is a quark and onto anti-quark flavors when a is an anti-quark. 
(For an anti-quark b, b denotes quark flavor.) The integration region in the convolution 
formula is specified by = and z^in = — d'ab-*c etc represents partonic 

hard cross sections which are the functions of Mandelstam variables in the parton level 
s = {xp + x'p')'^ = xx'S, i = {xp — Ph/zY = |T and u = {x'p' — Ph/zY = and are 
given as follows: 

-(iV^g + 2t)(g^ + -{NH + u-s)s 

O'ab^c = 0, 

- -Oab H TTTT^ OabOac, (H) 



ab—tc 
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^--^ = ^ ^4^^ ^ + ^^(iV^(5^ + 3S^n-2n^) + s3-s^n)5.„ 



0", 













yst 




4( 



+ \-A- + -]+-r{^^2 I, + ^U(12) 



Oan~^c. = < —Z —^^T-Z \ —TT—: > On 



^ N{m - l)Piu ' 

_/V2 

aan^n = ^ — f4s^ + lls^t + 19s^? + 22s^P + mH'^ + llst^ + 

+ N^^m-l)si^u^ {-St^^ + N\s^ + sH + 2sH- + si^ + t-)} , (14) 

where iV = 3 is the number of colors. A remarkable feature of (|T0|) is that the partonic 
hard cross sections for G'i7'(0, x) and Gpi'^iX) are the same, even though each diagram in 
Figs. 3~8 gives different contributions for the two functions. Accordingly, they appear 
in the combination of Gir(0,x) + Gi?(0,a;) in ffTOj) . Although we expect there is a reason 
for this simple result, we do not understand the origin at this point. Another important 
feature is that some terms in the hard cross section (the first terms in aab^c, ^al^ci (^ag^c 
etc.) have the same behavior as the SGP hard cross section which gives dominant 
contribution to An in the forward region {xp ^ I)- In addition their color factors are 
typically A^^ = 9 times larger than those in the SGP cross section, making the magnitude 
of the SFP hard cross section huge. This suggests that the SFP contribution gives rise to a 
significant contribution to the asymmetry at large xp even though the SFP function does 
not appear with the derivative. 

To see the impact of the SFP contribution, we have performed a numerical calculation of 
An for p^p nX, assuming that the SFP function is of the same order of magnitude as the 
SGP function. Kouvaris et al. [31] have parametrized the SGP function G^(x, x) so that 
the SGP contribution reproduces the A^ data obtained from RHIC and FNAL data. Their 
analysis shows that both data are reasonably well reproduced. In particular, they found 
that the derivative term in ^ brings dominant contribution compared to the nonderivative 
contribution. Here we adopt Fit(I) of [31] and assume G%{Q,x) + G^(0,x) = G%{x,x) 
(a = u,d). Fig. 9 shows A^ for the pion at \/S = 200 GeV and P^t = 1-5 GeV with 
and without SFP contribution. As is seen from the figure that the SFP contribution brings 
large effect in the positive xp region, while its effect is negligible in the negative xp region. 
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Figure 9: for p^p -> ttX at = 200 GeV and Phx = 1-5 GeV. Solid, fong-dashed, 
and dash-dot lines are, respectively, Aj^f for vr"^, vr" and 7r° obtained with only the SGP 
contribution. Dotted, short-dashed, and dash-double-dot lines are, respectively, Ajy for vr"*", 
7r~ and obtained with both SGP and SFP contributions. 

This is due to the large color factors for the terms in dab^c, o"ab^c) (^ag^c as noted 
above. From this figure it is clear that the SFP contribution can affect Ajy significantly 
even though it does not receive enhancement by the derivative unlike SGP contribution, 
unless the SFP function itself is small. The detail of the present calculation together with 
a phenomenological analysis of the data will be presented elsewhere |41j . 

To summarize, we have calculated the SFP contribution to the cross section for p^p -^■ 
hX associated with the twist-3 quark-gluon correlation functions in the polarized nucleon. 
We have also shown that its effect is significant and should be included in the analysis oi A^. 
For a more complete analysis, one needs to include the contribution from the triple-gluon 
twist-3 distribution function as well. Recent study [36] shows that for the open charm 
production has a potential to determine the function. We hope that the global analysis of 
all those data including all the QCD effects will clarify the origin of observed SSA. 
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